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I Abstract 

. We study Landau-gauge Yang-Mills theory by means of a nonperturbative vertex 

^ I expansion of the quantum effective action. Using an exact renormalization group 

equation, we compute the fully dressed gluon and ghost propagators to lowest non- 
. trivial order in the vertex expansion. In the mid-momentum regime, p"^ ~ 0(l)GeV^, 

OO I we probe the propagator flow with various ansdtze for the three- and four-point cor- 

' relations. We analyze the potential of these truncation schemes to generate a non- 

■ perturbative scale. We find universal infrared behavior of the propagators, if the 

, gluon dressing function has developed a mass-like structure at mid-momentum. The 

resulting power laws in the infrared support the Kugo-Ojima confinement scenario. 

1 Introduction 

^ , A remarkable property of four-dimensional Yang-Mills theories is the nonperturbative gen- 

^ I eration of a scale, separating the physics at high energies from that of low energies. Many 

■ - - ' high-energy phenomena can well be controlled with the elaborate machinery of perturbation 
theory. By contrast, essential low-energy properties such as confinement and bound-state 
formation as well as the transition region between the nonperturbative and the perturbative 
regime, are not yet fully understood. 

Genuinely nonperturbative frameworks such as functional methods for computing Green's 
functions or lattice Monte-Carlo simulations are required to investigate these phenomena. 
As a main advantage, functional continuum methods can cover orders of magnitude in 
momentum space and therefore naturally connect the high- and low-momentum regime. 
Analytical solutions are available in both the ultraviolet (UV) and also the infrared (IR) 
region of momentum. Certainly, the drawback of functional approaches to the infinite set 
of Green's functions is that one has to rely on truncations to obtain a closed, solvable 
system of equations. 

A prominent representative of the functional approach is the tower of Dyson-Schwinger 
equations (DSE) that has found a variety of applications in the context of QCD P1I21E1- In 
Landau-gauge Yang-Mills theory, the lowest order n-point functions, the ghost and gluon 
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propagators, have been investigated in various approximation and truncation schemes in 
the DSE approach [H HJ El 13 13 El • The picture emerging in the IR is that the 'geometric' 
degree of freedom, i.e., the Faddeev- Popov determinant, dominates over the dynamics of 
the gluon field: the ghost and gluon propagators are described by simple power laws with 
an IR finite or even vanishing gluon propagator (depending on the truncation scheme) and 
a ghost propagator more singular than a simple pole. In the momentum region accessible 
to lattice simulations, such a picture is nicely confirmed [TUl CH (121 UHl UH Uil ■ 

These IR power laws are in agreement with both Zwanziger's horizon condition arising 
from the gauge- fixing procedure and the Kugo-Ojima confinement criterion. The horizon 
condition is formulated as a boundary condition on the ghost and gluon propagators. It 
restricts the integration of gauge-field configurations in the Schwinger functional to the 
Gribov region, defined by a positive semidefinite Faddeev- Popov operator J^l- It has been 
shown that this restriction is sufficient to ensure that physical expectation values are not 
affected by gauge copies |17j. Entropy arguments have been employed to reason that the 
IR modes of the gauge and ghost fields are close to the Gribov horizon. The resulting 
horizon conditions state that the gluon propagator should vanish in the IR and the ghost 
propagator should be more singular than a simple pole [THj. 

The Kugo-Ojima confinement criterion (201 is a dynamical condition on the two- 
point function {D^cDyc), which ensures the conservation of global color charge. Provided 
BRST symmetry holds also nonperturbatively, one can then define a physical state space 
Vphys containing colorless states only. In Landau gauge the Kugo-Ojima criterion can be 
translated into a condition for the IR behavior of the ghost propagator: it is fulfilled if the 
ghost propagator is more singular than a simple pole |i2j. What remains to be shown in 
this scenario is the appearance of a mass gap and the violation of cluster decomposition in 

Functional methods for computing Green's functions can be implemented in various 
different but interrelated formulations, each with its own assets and drawbacks. In this 
work, we employ the exact renormalization group (RG) [221 12H1 121] formulated in terms 
of flow equations for the n-point functions. The exact RG flow is derived from the scaling 
behavior of the effective action with respect to an IR cutoff k. Starting with an initial 
condition for the effective action at a UV cutoff scale Auv, e.g., in terms of the bare action, 
the RG flow equation is integrated from k = Auv down to k = 0, thereby generating the 
full quantum effective action. If the effective action is expanded in powers of fields, these 
couplings are the (IPI) Green's functions (proper vertices) of the theory. 

The RG equations are, for instance, intimately related to the DSEs.^ On a formal 
level, it can be shown that effective actions satisfying the DSEs in the presence of an IR 
cutoff k are (quasi-)fixed points of the exact RG flow. Thus a solution of the flow equation 
for — > is also a solution of the ordinary DSEs ITTj . This property should hold 
in an approximate sense even if the effective action is truncated, provided the truncation 
preserves the dominant structure of the theory. 

^There is also a close relation to the functional approach based on 2PI (or, more generally, nPI) effective 
actions which are currently investigated in the context of non-equilibrium quantum field theory [21] . 
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In comparison to the corresponding DSEs, the flow equation for the Yang-Mills prop- 
agators has been much less explored. The first investigation was reported in fIE\ 1^ . 
From the behavior of the propagator functions around Aqcd, a strong IR divergence of the 
gluon propagator consistent with the old idea of infrared slavery was conjectured. Similar 
observations were made in [22] • These results seem to be in marked contrast to the IR 
power laws described above; however, in these works, the integration of the flow was only 
performed down to a finite value of k, and the deep IR — > was not explored. Indeed, a 
recently developed flow-equation technique for performing a self-consistency study of the 
IR scaling of propagators flnds power laws in agreement with the ones obtained in the DSE 
formulation, thereby reconciling the two approaches 

In this work, we employ the RG flow equation in order to study the full flow from the 
perturbative UV to the deep IR, investigating the possible mechanisms that connect the 
high-energy behavior of the propagators with their IR power laws. In particular, we will 
not use self-consistency criteria for determining the IR power laws, but solve the flow in 
order to see if and how the power laws emerge. The paper is organized as follows: in the 
following, we flrst summarize aspects of Landau-gauge Yang-Mills theory that are relevant 
to our work. In section El we introduce the exact RG equations for the propagators of 
Yang-Mills theory and discuss our initial truncation of the effective action. The role of 
gauge symmetry in controlling the flow is specifled. In section HI we flrst discuss important 
qualitative features of the flow equation and then present our numerical results. We focus 
in particular on the delicate mid-momentum regime, p"^ ~ C(l)GeV^ by probing the 
propagator flow with various vertex ansdtze. Finally, we solve the flow towards the IR, 
in order to explore the "domain of attractivity" for which the IR power laws represent an 
IR stable flxed point. This again provides information about the dynamical mechanisms 
required in the mid-momentum regime, which still remains insufficiently understood. We 
end with a detailed discussion of our results in the conclusions. 

2 Propagators in Landau-gauge Yang-Mills theory 

In the continuum Green's function approach to Yang-Mills theory. Landau gauge has been 
a favorite choice for a number of reasons. First, of all linear covariant gauges. Landau 
gauge is the only one symmetric under the exchange of ghosts and antighosts. This, on 
the one hand, allows one to interpret ghosts and antighosts as (unphysical) particles and 
antiparticles. On the other hand, this symmetry is a useful guiding principle for the 
construction of a nonperturbative ansatz for the ghost-gluon vertex [HI E]- As will be 
detailed later on, such an ansatz is necessary in both the DSE and the exact RG flow- 
equation framework in order to close the associated equations determining the ghost and 
gluon propagators. Secondly, the ghost-gluon vertex does not suffer from UV divergencies 
in Landau gauge |S2]; i-e., the vertex renormalization factor Zqj^^q can be set to one. Again, 
this provides a useful constraint on possible nonperturbative vertex dressings. In fact, as 
has been shown in the DSE framework ^51 El 13 El , even the bare ghost-gluon vertex is a 
good approximation in the UV and IR. 
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Thirdly, a direct consequence of UV finiteness of the ghost-gluon vertex in Landau 
gauge is a nonperturbative definition of the running couphng in terms of the propagators. 
In Euchdean momentum space, the ghost and gluon propagators are given by 

Dac(p-) = (1) 

with the ghost dressing function G{p'^) and the gluon dressing ZijP')? With the help of 
the STI, 

1 = ZcAC = ^y^5 (3) 

relating the vertex renormalization factor with the corresponding factors for the coupling 
g, the ghost and the gluon fields, we can define a running coupling a(p^) by [221 El 

a{p') = ^G'{p')Z{p'). (4) 

Note that the right-hand side of this definition is an RG invariant, i.e., a(p^) does not 
depend on the renormalization scale. Within the Green's functions approach, the IR 
behavior of the ghost and gluon dressing functions can be determined analytically. For 
momenta p -C Aqcd, the dressing functions are given by simple power laws, 

Gip') ~ (/)~^ Zip') ~ (5) 

with interrelated exponents [H El El 13 El EHj- Hereby k is an irrational number, n = 
(93 — vl201)/98 ~ 0.595, which depends slightly on the truncation scheme [Hj. From the 
power laws, Eq. (0), we see immediately that the running coupling, Eq. (HJ, has a fixed 
point in the IR. 

A fourth reason why Landau gauge is interesting has already been mentioned in the 
introduction: in this gauge, there is a direct connection between the Kugo-Ojima confine- 
ment scenario and the ghost propagator. If the ghost propagator is more singular than a 
simple pole, global gauge symmetry is unbroken and one can demonstrate that the physical 
state space of the theory, defined as the cohomology of the BRST operator, contains color 
singlets only ^T] . 

Finally, Landau gauge is known to be a fixed point of the renormalization flow jHlj. 
Within a given truncation, it is possible to show that this fixed point is IR attractive for 
a wide range of initial gauge parameters |2E1- This suggests that an investigation of Yang- 
Mills theory starting in a general linear covariant gauge at large IR cutoff will generically 
end up in the Landau gauge after the flow has been integrated. This, together with the 
simplifications mentioned above, speaks for Landau gauge as a natural and convenient 
choice from the very beginning. 

^Here we employ the conventions of the DSE hterature; in the RG hterature, the dressing functions are 
usually inversely defined: e.g., Z{p'^) 1/Z{p'^). 
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3 Flow equation for the vertex expansion 



3.1 Exact renormalization group 

To provide a short summary of the exact RG approach, let us begin with the gauge-fixed 
action for SU(A^c) Yang-Mills theory in a covariant gauge in D dimensional Euclidean 
spacetime, 



"S'ym — / dPx 



(6) 



In the fiow equation approach, an IR cutoff k for the quantum fluctuations is implemented 
by adding a GauBian term to the action. 



ASk = - I d'^xAlR'l^^Al + / d^xCR^^CK (7) 



The momentum-dependent regulators Ra,c cut off the IR fluctuations for gluons and 
ghosts, respectively, at a momentum scale ~ k"^. Defining the quantum theory via 
the Schwinger functional with the action S'ym + AS'^, a Legendre transform leads us to the 
effective average action r^. It already contains the effects of all quantum fluctuations with 
momenta larger than fc, and governs the dynamics of the remaining modes with momenta 
smaller than k. The response of the effective action under a variation of the cutoff scale k 
is described by the flow equation, 

dtTk[<p\ = ^Tr6?A[0]9ti?A - TiGcWtRc. (8) 



where (f) = {A, C, C) and t = In k/A. The traces run over all indices including momentum, 
and Ga/cI'P] — i^k M+-^)a/c abbreviates the full regularized gauge/ghost fleld propagator 
with r^.^"* [(f)] = 6^Tk[(j)\/S(j)S(j). Once initial conditions are specifled at a high scale A in terms 
of the microscopic action to be quantized, Tk=A — * S'ym, bare, the flow equation describes 
the RG trajectory of the effective average action towards the IR. The end point F = Tk=o 
corresponds to the full quantum effective action. The flow equation is UV and IR flnite 
by construction, it is exact despite its one-loop structure, and the momentum trace is 
dominated by momenta ~ fc^. 

For a good approximate solution of the flow equation, the ansatz for the effective 
action should include the relevant degrees of freedom for the problem under consideration. 
In this work, we assume that an expansion of the effective action in terms of full vertices 
can illuminate aspects of the nonperturbative structure of the theory. 



^ n 



S{P1 + ■ ■ ■ + Pn) rl^^Pl, ■■■,Pn) <P{Pl) ■ ■ ■ (PiPn), (9) 



where = / d^p/{2n)^, and 6{p) = {2t:)^ 5^^\p) . Inserting Eq. © into Eq. © and 
taking appropriate derivatives, we obtain an inflnite set of coupled flrst-oder differential 
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equations for the proper vertices F^"''. Truncating the expansion at order nmax leaves all 
equations for the vertices r^^"-"™'''' ^) unaffected. In order to close this tower of equations, 
the vertices of order nmax and rimax — 1 can either be derived from their truncated equa- 
tions or taken as bare - or even built upon inspired guesswork. This defines a consistent 
approximation scheme that can in principle be iterated to arbitrarily high orders in rimax- 

3.2 Truncated flow equations 

In this work, we set n^ax = 4 and solve the dynamic equations for the 2-point vertices: 
the inverse propagators. To this end, we employ either bare 3- and 4-point vertices or use 
constructions derived from further physical information. The tree-level- vertex structure of 
the Landau gauge is exploited, e.g., by setting the four-ghost vertex to zero. Our truncation 
for the effective action is given by 



q 



+m I + q2 + qs) VS% uili^ q2, qs) C^i-qM'MC'iqs) 



|ji>92,g3 



+i^?A / Siqi + q2 + qs) V.t^.iqu q2, qs) Al{q^)Al{q,)Al{q,) 

+-,91 1 Hqi + ■■■ + q^)v,Tiqu • • • , qMU^Mliq2)A''^{qs)Atiq4), (lo) 

^ Jqi,...q4 

where qa denotes the renormalized gauge coupling at the high scale A where our flow will 
be initiated. In the first line, the fc-dependent dressing functions Z{p'^), Zi^^p^), GijP') char- 
acterize the fluctuation-induced modifications of the bare propagators of the transversal 
and longitudinal gluons and the ghosts, respectively. These are at the center of interest in 
the present work. Furthermore, a gluon mass m\ has been included explicitly, the role of 
which will be explained in detail later on. 

Inserting this truncation into Eq. (jHl), the flow equation reduces to a set of equations 
for the propagators which we display diagrammatically in Fig. ^ In the following, we will 
consider truncations where the Lorentz and color index structures of the vertices are kept 
identical to the bare ones: 

K?AC,M(^i'^2,g3) = r"'qi,VcAc{quq2.qs) 

Vs%i<lu<l2,q3) = r'''qi,V^A{qi,q2,qs) 
V,T{qi,...,q,) = r'r'V,A{qu...,q,). (11) 

Here V^^c^ ^sa, are functions of the momentum invariants of g^^, examples of which 
will be given later; the bare-vertex truncation is defined by Vc^c = ^sa = V^a = 1 up to 
RG rescalings of the fields (cf. Appendix IH]) . 

Even though exact solutions of the flow equation are independent of the choice of the 
regulator, approximate solutions do generally depend on the specific form of R. A careful 
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Figure 1: Flow equations for the propagators. Filled circles denote full propagators, open 
circles correspond to full vertices, and crossed circles indicate the (9t-Rfc-operator insertion. 

choice of the regulator improves numerical stability as well as the quality of the truncation 
in general Here we concentrate on regulators of the type 

Ra{p') = (^^t(p) + ^Aip'/k'), Rcip') = ^^rcipVe). (12) 

The dimensionless regulator shape function rA,c'{x) with a dimensionless argument de- 
termines the analytic form of the momentum space regularization. In addition, we have 
included the dressing functions in the regulator. In the first place, this construction guar- 
antees that the flow equations are invariant under RG rescaling [30]. Moreover, since we 
expect a fast RG running of many couplings in the nonperturbative domain, the inclusion 
of the dressing functions accounts for an adjustment of the regulator to the spectral flow 
of quantum fluctuations [201 IHIl IHHI • 

Let us now introduce some notation in order to facilitate a more explicit representation 
of the flow equation for the propagators. First we specify squared momentum variables, 
X, y, z, involving the external momentum p and the loop momentum q, 

x = p^, y = q^, z = {p- qf = x + y -2u^, {x,y, z) = {x,y, z)/k'^, (13) 

where u := p.5'/(|p| |5'|) denotes an angle variable. Here and in the following all hatted 
quantities are dimensionless by appropriate rescaling with k^. Although we are interested 
in D = 4 dimensional spacetime in this work, it is worthwhile for future application to 
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perform the computation in arbitrary D dimensional spacetime. For this, we need the D 
dimensional loop integral of an arbitrary function of x, y, u: 



^ f{x,y,u) = vd [ dyy° ^ [ du{l-u^)°^'' f{x,y,u), 



(27r) 



D 



(14) 



^0 J -I 

2 r(^) 1 



Here f d is related to the surface volume of a -D dimensional sphere ^. In order to display 
all dimensionful quantities in terms of fc^, it is also useful to introduce the dimensionless 
coupling ^A, 

9l = 9lk''-\ (16) 

where g\ denotes the value of the coupling at an initial high scale A. 

After insertion of this specific regulator form given above, the flow equation for the 
inverse transverse gluon propagator r^^^^(p^) = / Z{p'^) + m\ can be written as 



dt 



X 



Zix] 



+ m] 



(*00 pi 

-AglN.vok^ j dyy^-^ J du{l-u^' 

Q^'^ix, y, u) K^lfix, y, u) V^a Z{ye)Z{ze) V^a 

+g«'^(x, u) Kgfix, y, u) VcAc G{yk^)G{ze)VcAc 
-^^^^glKvn e dy yT-i (17) 

( " ^'""^^^^ " \yTA{y)d^nZ{yk')^ V,A. 

where the last term is the tadpole contribution. Since we are ultimatly aiming at the 
Landau gauge, we have omitted any contribution of the longitudinal modes on the right- 
hand side of the flow equation which decouple in this limit. The flow equation for the 
inverse ghost propagator r^^^p(p^) = /G{p'^) then yields 



g-3 
2 



a 



X 



G{x) 



dyy^-'J , 



du [1 — 



Q^'^ix,y,u)K^,fix,y,u) 



xVcAcG{ye)Z{ze)VcAc 



:i8) 



The quantities Q^'^, Q'^''^ denote the kernels of the gluon and ghost loop, respectively, in 
the gluon flow equation, whereas Q*^'^ represents the ghost equation kernel. These kernels 
are given in Appendix ^ 
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The quantity K^eg abbreviates the contribution from the regularized propagators which 
reads: 



'reg ^{^iVi'^) 



1 



Pa{y)Pb{z) 



-y'^Kiy) - lyraiy)dt\nZaiyk^) 



+ 



Pb{z) 



(19) 



where the prime denotes a derivative with respect to the argument. Here Za^b denote gluon 
and ghost dressing functions, Z^p^) and G{jP'), and pa,b abbreviate the corresponding 
regularized inverse propagators, 



Pc{y) = y (1 + rc{y)), PA{y) = y (l + rA{y)) + Ziyk"^) ml/k^. 



(20) 



An explicit expression for the regulator shape function can be found in Appendix iBl 

Note that all scale derivatives dt = k{d/dk) that occur are taken at fixed dimensionful 
external momentum, dtx = = dtlxk"^), so that the scale derivatives on both sides of the 
flow equation have the same meaning. However, on internal momentum variables the scale 
derivative acts on the manifest fc^ dependence, dty = 0. This implies the somewhat clumsy 
differentiation rule for the dependent momentum variable dt{zk'^) = 2yk'^ — 2u^yyk'^ 



Equations (fT7|) and (fTH|) represent a closed set of equations for the gluon and ghost 
propagators in the Landau gauge. Once we have specified initial conditions for Z and G 
at, say, a perturbative UV scale = A, we can integrate the flow equations down to A; = 
and read off the form of the fully dressed propagators. However, before we do so, we first 
have to discuss the issue of gauge invariance. 



3.3 Gauge invariance 

Since we are working in a gauge-fixed formulation of Yang-Mills theory, gauge invariance 
of the system is encoded in a constraint for the effective action. This constraint can 
either be formulated as a Ward-Takahashi identity or, invoking the BRST formalism, as a 
Slavnov- Taylor identity. In addition to the breaking of gauge invariance by the gauge-fixing 
procedure, the regulator term ((7j) represents another source of gauge or BRST symmetry 
breaking. In order to account for this additional breaking, both the Slavnov- Taylor and the 
Ward-Takahashi identity are modified by further terms depending on the regulator. These 
guarantee the restoration of gauge invariance in the limit k — *• 0.^ The BRST formalism in 
the flow equation framework has been worked out in Refs. [211 EHl 1^ ■ Since the resulting 
equation is important for the results of the present paper, we rederive the same findings 
from the Ward-Takahashi identity here, following 

Using the momentum-space representation of the generator Q°- of infinitesimal gauge 
transformations formulated in terms of unrenormalized coupling and the fields Aq^, Cq , Cq 

^ As an alternative to the present formalism, the construction of manifestly gauge-invariant flows has 
been put forward in |39|. 
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which occur in the microscopic action S, 

6 



(21) 

the modified Ward-Takahashi identity (mWTI) for covariant gauges can be written as"^ 



Q'ip) 



T9or'' I iP + Q),iP + qUAl^,{-q)Al,ip + q))c 



q 



-9or'' I pMp + ?i - 92)92^^"^ - -mf'^Al^ip + q,- q^)) {Qiqi)CRq2))con 

>^ 91:92 

_ig2fatcffet I p^l^Q^q^^Al^^ip + gi - q2)C^{q2))con 

>^ 91:92 

-gor"' [ [Rcip + q)- Rc{q)] {C^oi^Kiq + p)),on, (22) 

Jq 

where the first four fines represent the standard Ward-Takahashi identity, and the last 
two fines are tfie modification owing to the regulator. In the limit — *• 0, the regulator 
terms vanish such that the standard WTI is recovered. The important observation is 
that the mWTI is compatible with the fiow equation in the sense that if a solution of 
the fiow equation satisfies the mWTI at one scale, it does so for all scales. Hence, if our 
initial conditions at the high scale satisfy the mWTI aX k = A, the end point of the RG 
trajectory at A; — > which is the full quantum effective action will be gauge invariant, i.e., 
will fulfill the standard WTI. However, these statements only hold for the full effective 
action. A truncated flow can leave the gauge-invariant trajectory which is constrained by 
the mWTI. In order to satisfy the symmetry principle encoded in the mWTI, the degrees of 
freedom in the truncation can be subdivided into truly independent ones, the dynamics of 
which is determined by the flow equation, and the dependent ones that can be expressed in 
terms of independent degrees of freedom by virtue of the mWTI. In this way, the symmetry 
principle is satisfled on the theory subspace deflned by the truncation, and gauge invariance 
is consistently implemented within the truncation [101 H^ . 

In the present truncation, the mWTI flrst of all constrains the longitudinal gluon prop- 
agator. Although the longitudinal components decouple in the Landau gauge which we are 

*Here we use a convenient representation of the mWTI in terms of connected n-point functions to be 
evaluated from the Schwinger functional in presence of the regulator term 0; of course, by Legendre 
transformation, the connected n-point functions can be expressed in terms of the effective action Tk[(j)] 
and derivatives thereof. Moreover, our notation does not distinguish between the fluctuation fields (to 
be integrated over in the Schwinger functional) and the corresponding "classical" field (conjugate to the 
source), but the meaning should be obvious from the context. 



10 



aiming at, the mWTI nevertheless provides an important piece of information. Evaluating 
Eq. ()22j) in the zero-momentum limit of the longitudinal propagator, the mWTI constrains 
the value of the gluon mass m^. In order to represent this constraint for the gluon mass 
in terms of the quantities appearing in our truncation (fTIUl . we need to know the relation 
between the fields {Ao,Co,Co) appearing in Q"- and the fields renormalized at the scale A 
in which the effective action is expressed, i.e., we need to determine the renormalization 
factors Za and Zc defined by 

Aof. = ^A^, {Co,Co) = ^/Z^{C,C). (23) 

The connection between the renormalization factors in the exact RG approach and the 
conventional MS scheme have been worked out in Ref. jH]. Identifying A with the renor- 
malization scale, the boundary conditions for the flow at A are equivalent to renormalization 
conditions in the MS scheme, and the finite parts of the respective renormalization factors 
can be related to each other. 

As will be detailed in Sect. 14.11 our initial effective action at A will be obtained by 
integrating the bare action Syu + AS^, Eqs. (j6l7p . from an even larger scale Auv down to 
A. With 

1 = Zr,A,v(0) = ^a(0)^a(A), 1 = Ga,,(0) = Ga(0)Zc(A), (24) 
we then obtain 

^^ = 7^1(0)^0^, iC,C) = y^G^iCo,Co), (25) 

where (Za(p^), G'a(]9^)) = \imk^A{Z{p'^),G{p'^)). Furthermore, the nonperturbative defini- 
tion of the running coupling, Eq. (^, gives the relation 

9'o=9lGliO)ZA{0). (26) 
With these relations the mWTI for the gluon mass yields 

< = ^^-om Ga(0) e r dy y-''-'U,c G(yk^)^-=W^ + (D-l/'' 



X 




fr'^iy) Z{yP)ml{fr'A{y) +yrAiy)) , Z'{yP)mlfrAiy) 
yrA[y) + —FT^ u Z-77~\ ^ 



4 paH)) F pA{y) A;2 p^(y) 

(27) 

where we have again dropped longitudinal contributions on the right-hand side according to 
the Landau gauge, and primes denote a derivative with respect to y. The occurrence of the 
gluon mass ruk is an obvious consequence of gauge-symmetry breaking by the regulator. 
Assuming that the dressing functions diverge, if at all, with a simple power or weaker 
V-iAZip^), VcAc^ip^) ^ (l/p^)" in the IR limit, the k dependence of the mass is bounded 
by m^, < {k'^Y~'^ towards the IR. Removing the regulator with A; — ^ 0, the mass vanishes 
as long as K < 1. Therefore, k < 1 serves as a consistency condition that our solutions 
have to obey in order to satisfy the standard WTI in the present truncation. 
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Since the gluon mass also appears on the right-hand side of Eq. (P7j). the mWTI rep- 
resents a "gap" equation for a self-consistent determination of the gluon mass. Moreover, 
the RG evolution of this mass is completely determined by the evolution of Z and G and 
the vertices. In this way, the mWTI guarantees that the gluon mass is not an RG relevant 
parameter, as is conventionally the case for boson masses, but RG irrelevant, since it drops 
to zero with ^ 0. 

On the other hand, from the flow equation for the transversal gluon propagator ()17|) . we 
can read off the flow of the gluon mass in the zero-momentum limit. Of course, for the exact 
solution, this mass flow would fully agree with that of the mWTI; both mass equations 
also agree in the perturbative limit of the truncation (G, Z const., Voacaam ~^ 1); 
which serves as a nontrivial check of the approach j^EI- Nonperturbatively, however, both 
equations for the mass can differ strongly, even qualitatively, since an exact implementation 
of the gauge constraint always requires information which is beyond the truncation. The 
gluon mass is particularly sensitive to this fact, since any deviation from a gauge-invariant 
RG trajectory can lift its protection against quadratic renormalization (as can occur for 
massive bosons); then, the gluon mass would be of the order of the cutoff, ~ A^. Indeed, 
this is what would happen if we took the mass flow of the transversal gluon seriously, 
disregarding the mWTI. In similar truncations of DSEs, the very same fact occurs as 
quadratic divergencies in the gluon equation which have to be removed by hand. In the 
context of flow equations, we have the mWTI at our disposal for controlling the mass flow. 
For this, we subtract the mass flow from the transversal gluon equation and put in the 
mass from the mWTI: 

= RHS C3L.(^wTi) -^^^'Ifiow = C3L,(^^^j^ -RHS CB|_o, (28) 

where the gluon mass in the RHS is taken from the mWTI. For the exact solution, Eq. ()28|1 
corresponds to a zero operation, whereas in the truncation, this procedure removes quadrat- 
ically renormalizing parameters by virtue of the gauge constraint.^ 

4 Flow analysis 

In this section, we analyze the flow equation partly analytically but mainly numerically for 
various initial conditions and truncations. We concentrate on the case of D = 4 dimensional 
spacetime and Nc = 'i colors unless stated otherwise. In the next subsection we will specify 
two different types of boundary conditions at large and intermediate momenta that serve as 
initial conditions for the flow. We proceed by analyzing the RG equations analytically on a 
qualitative level and determine conditions for the higher n-point functions such that the IR 

^In |2ni, the difference of the propagator flow resulting from using either the mWTI mass flow or 
the direct mass flow was taken as a dynamical criterion for the breakdown of the truncation. In view 
of the qualitatively different RG behavior of the two masses in this truncation, this criterion is rather 
restrictive. The flow mass corresponds to an unphysical quadratic divergence in this truncation and hence 
can be expected to deviate signiflcantly from the mWTI mass. Here we abandon this criterion and give 
preference to the mWTI mass. 
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attractive power-law domain can be approached by the flow. Guided by these findings, we 
then perform a numerical analysis of the fiow in various truncation schemes. The details 
of our numerical methods are given in Appendix IdI 

4.1 Initial conditions 

The solution of a quantum theory is determined by the RG fiow equation and the cor- 
responding initial conditions. For the latter, we have to specify the details of the initial 
effective action Fa, serving as our boundary condition at the scale A. For a perturbatively 
small coupling g;^ at A, one could choose the bare action (jH]) for Fa. Improved initial 
conditions can be obtained by perturbatively integrating the flow equation from a larger 
scale Auv down to A, employing bare dressing functions on the right-hand side of the flow 
equation. This guarantees the correct one-loop behavior of the dressing functions and sub- 
sequently of the running coupling for momenta ^ A far larger than our starting scale. 
Of course, a correct matching to the full flow at A requires the inclusion of the regula- 
tor function in this perturbative determination of the initial conditions. Qualitatively, 
the high- momentum behavior, ^ A^, of the dressing functions for one- loop improved 
perturbative initial conditions (PIC) is 



where 7a = —13/22 and 7c = —9/44 denote the anomalous dimensions for gluons and 
ghosts, respectively, for arbitrary A^c, and Z\,Ga = Z(A^), G(A^)|a denote the normaliza- 
tions of the flelds. For low momenta ^ A^, the perturbative initial conditions for the 
dressing functions go to a constant, owing to the presence of the regulator.^ 

Certainly, universality tells us that such a one-loop RG improvement is irrelevant for 
the questions related to the IR behavior. Whether a truncation is good enough to stabilize 
the flow in the mid-momentum regime between the perturbative UV and the deep IR 
should not depend on the details of the UV behavior of the initial effective action. We 
have conflrmed that this is indeed not the case, employing both a bare initial action and 
its one- loop improved counterpart. 

Since it turns out that a full solution connecting perturbative initial conditions with the 
expected IR behavior is difficult to construct, we design instead a set of mid-scale initial 
conditions for exploring in particular the fiow in the nonperturbative IR regime. For this, 
we consider the initial scale A to be already in the nonperturbative regime but still high 
enough such that the IR asymptotics has not yet developed, say, for instance, A > IGeV. 
For definiteness, we study two contrary cases: one in which the ghost dressing function has 

^In principle, we could also use improved initial conditions in the form of one-loop leading-log resummed 
perturbation theory. However, such an improvement would hardly exert any influence on the remaining 
flow towards the IR, which is the focus of the present investigation. 




(29) 
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already developed a scale L, whereas the gluon remains perturbative (or even constant), 



) 



vc 



A 



Z(''^)(p2)|^ = ZA. 



(30) 



Here, rjc denotes a trial exponent to be varied, and the scale L separates perturbative from 
deep IR behavior. The other case is given by a scale in the gluon dressing function and a 
perturbative ghost. 



For a trial gluon power 77^ ~ 1, the scale L plays the role of a nonperturbative gluon mass. 
Such a nonperturbative gluon mass should not be confused with the mass controlled 
by the WTI, cf. subsection 13.31 In general, the scale L will survive when the IR cutoff k 
is taken to zero, whereas ^ in this limit according to the WTI. 

Note that Eqs. (jHUj) and (jHT|) are initial conditions at a fixed mid-momentum scale, and 
do not represent an ansatz for the form of the dressings at even lower scales. In other 
words, we will determine the resulting full momentum dependence at lower scales and not 
merely the flow of the parameters rjA^c and L. 

The use of these mid-scale initial conditions (j3Up . (j31|) for the IR flow facilitates a 
study of the IR attractive domain of the asymptotic power solutions obtained so far in 
the literature. In this way, we can analyze the necessary ingredients for a solution of the 
flow at an intermediate scale k to finally run into the power-law solutions for A; — 0; it 
dispenses us from the need to solve the complete flow from the UV to the IR in particular 
in the transition region, where vertex corrections are expected to become important. 

4.2 Qualitative analysis 

Some important features of the present set of RG equations for the propagators can already 
be extracted by a qualitative analysis. We begin with the observation that the perturbative 
initial conditions (j29|) are such that the dressing functions are increasing functions for 
decreasing momentum = x. In view of the expected IR behavior G{x) ~ x~'^ with 
K > 0.5 (cf. Sect. El), this property should persist and even be enhanced for the ghost 
dressing. By contrast, the gluon dressing is expected to decrease in the IR with Z{x) ~ x'^'^, 
such that the flow equation has to bend the gluon dressing in the IR. 

Let us check the ghost equation (|18|) first. As long as the vertex dressing Vq^^q remains 
positive, the RHS of the ghost flow is strictly positive; particularly, the kernel Q'^'^, given 
in Appendix ^ and the dressing functions are positive. The first term in the regulator 
quantity Kj-eg is positive as r'{y) < for all monotonous regulator shape functions, like 
the one specified in Eq. ()D.7|) below. Moreover, the mechanism for a sign change should 
not depend on the details of the regulator. The other terms ~ dt In Za in K^^g are positive 
for perturbative initial conditions, and thus they cannot induce a sign change, owing to 
the dt structure. Regarding the RHS of the ghost equation as "/? functions" for a set of 




A, 




(31) 
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"couplings" x/G{x) labeled by x, the positivity of these "beta functions" guarantees that 
the "couplings" decrease towards the IR. Consequently, this decrease of x/G{x) with k 
implies an increase of G{x) towards the IR. This is consistent with the IR expectation of 
an enhanced ghost, G{x) ~ x~'^. 

Next we apply this argument to the gluon equation (fTTjl . Since we anticipate a suppres- 
sion of Z{x) for small x, we need a negative "/3 function" for x/Z{x). Hence, we expect 
a negative RHS for small x, once we have subtracted the gluon mass flow dtml in order 
to isolate the flow of x/Z{x). Let us analyze Eq. (fT7j) term by term, beginning with the 
gluon loop with kernel Q^'^ . For simplicity, we consider first the limit of x <^ where 
the propagator and vertex dressings as well as i^^^eg become rather independent of x, owing 
to the regularization (the z dependence can safely be replaced by y in these quantities). 
Hence the sign of the kernel after mass flow subtraction determines the overall sign. Using 
representation ()A.1|) . this corresponds to 



Q^'^ix, y, u) - Q^'^(0, u) = -(1 - u') 



X + 



xy 



1 



u 



D-1 



< 0. 



(32) 



Together with the overall minus sign in Eq. ()17|) . we find that the gluon loop contribution 
to the "/5 function" is positive. Hence x/Z(x)|giuon loop decreases towards the IR, implying 
that ^(x)lgiuon loop increases. We have derived this result in the limit x <^ k"^; but owing 
to the regularization, a well-converging expansion of the flow in x/k"^ exists, such that the 
validity of this statement can be extended to x < fc^. Even within the limits of these mild 
assumptions, we arrive at an important first result: the gluon loop in the gluon equation 
cannot be the source of a bending and a subsequent suppression of the gluon dressing. 
This statement holds for all 3-gluon vertex dressings that are non-negative and preserve 
the bare-vertex index structure. 

Now we turn to the ghost-loop contribution to the gluon equation with kernel Q'^'^ . 
The situation here is more subtle, since the kernel is independent of x (see Eq. ()A.2|) ). 
such that the residual terms after mass-flow subtraction arise from the i-dependent ghost 
and vertex dressings and -ft"reg- Collecting all these z dependencies in a function f{z) ~ 



V(j j^(jG{zk )K^c^, we have to study the sign of 



f{z) - f{y) = xf'iy) + 0{x' 



(33) 



where we have dropped all terms odd in u that are canceled by the u integration. G{p'^) as 
well as iiTj-eg for generic regulator shape functions are decreasing functions of the momentum. 
In the bare- vertex truncation Vj^^,^ = 1, the RHS of Eq. (j33|) is therefore negative for small 
X, because f'{y) < 0. Together with the overall minus sign of Eq. (fTTj). this again implies a 
positive "/? function" and thus a contribution to the flow of Z{x) that cannot suppress the 
dressing function in the IR. This statement is, of course, weaker than the one given above 
for the gluon loop, since a non-negative ghost-gluon vertex dressing can, in principle, turn 
the sign of f\y). However, this would be an unexpected source for the bending of the gluon 
dressing, since the ghost loop is not anticipated to play a dominant role at intermediate 
scales. 
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Let us finally continue our discussion with the tadpole contribution in Eq. (fTTj). For 
a bare- vertex truncation V^a const., the tadpole does not contribute at all to the flow 
of Z{x), but only to the mass flow. However, as soon as the 4- vertex acquires a nontriv- 
ial momentum dependence under the flow, the tadpole can contribute with either sign. 
Especially if the vertex dressing decreases with external momentum x, its "/5 function" 
contribution can be negative, thereby inducing a suppression of the gluon in the IR. 

To summarize, a vertex expansion of the effective action has the potential to describe 
Landau-gauge ghost enhancement and gluon suppression in agreement with, e.g., lattice 
results. A realization of this scenario in a truncation on the 3-vertex level requires excep- 
tional (i.e., negative) vertex dressings or suitable index structures beyond the bare ones. 
By contrast, a truncation including dressed 4- vertices can accommodate gluon suppression 
more conventionally with merely positive vertex dressings and bare index structures.^ 

4.3 Vertices 

Since we intend to perform a vertex expansion to lowest non-trivial order, we do not 
attempt to solve the full dynamical equations for the vertices. The choice of the latter 
hence determines the remaining unspecified part of the truncation. Basically three different 
strategies can be pursued: first, the highest vertices in a vertex expansion can be taken 
as bare. Second, these vertices can be computed self-consistently from their truncated 
dynamical equations. These first two strategies define a consistent approximation scheme 
generalizable to higher orders. Third, the vertices can be modeled employing further 
physical information. 

In the literature, the second strategy of deriving or at least constraining vertices by 
their truncated equations has been frequently followed. In particular in gauge theories, the 
constraints from the Ward-Takahashi or Slavnov- Taylor identities can serve as an additional 
input and have been exploited in fl^ |2Hl H] • However, it still remains unclear whether this 
strategy can be successful, since the truncated equations may lack important structural 
information from the neglected terms; in fact, this manner of construction can even lead 
to inconsistencies as exemplified in IHl EHl • 

In this work, we use a set of different vertices in order to explore different routes that 
the flow can take from the UV to the IR. First of all, the ghost-gluon vertex is taken to 
be bare, Vqaq = 1, which is in agreement with the non-renormalization property but 
neglects possible finite renormalizations. In the gluonic sector, we consider: 

(a) bare 3-gluon vertex, V^a = Ga/Za, where factors G\,Z\ take the normalization of 
the fields at A into account (cf. Appendix : 

(b) modified bare 3-gluon vertex, Vsa = pG\/Z\, where the possibly /c-dependent pa- 
rameter p can be varied in order to suppress or enhance the 3-gluon vertex; we will 

^Yet another option by which the observations of this section could be evaded is given by the freedom of 
performing fc-dependent field redefinitions under the flow |45j . This property of the flow can be exploited 
in order to optimize the degrees of freedom in a truncation. However, in the absence of a convincing 
optimization criterion, we will not explore this possibility any further in the present work. 
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especially make use of a vertex suppression with p becoming gradually smaller than 
1, in agreement with recent lattice results ^7j. An exceptional case is given by p = 0: 
the "ghost-loop-only" truncation; 



(c) dressed 3-gluon vertex, 



1 1 



(34) 



Vaaa 



where and {p — qY denote the momenta running around the loop, and 5 = —9/44 
is the anomalous dimension of the ghost dressing function to one-loop order. This 
type of vertex construction has been used in the DSEs to ensure the correct one-loop 
anomalous dimensions of the ghost and gluon dressing functions in the UV . Since 
it breaks RG scaling invariance explicitly in the nonperturbative momentum regime, 
rescaling the initial conditions for the fields with G\, Z\ corresponds to suppressing 
or enhancing this vertex analogous to (b); 

(d) bare 4-gluon vertex V4A = G\/Z\. However, since the 4-gluon vertex occurs only in 
the tadpole diagram, the bare-4-gluon vertex truncation is identical to the no-4-gluon- 
vertex truncation, since the momentum-independent part of the tadpole contributes 
only to the mass flow and thus is subtracted completely; 

(e) dressed 4-gluon vertex. 



where p"^ denotes the external momentum, and can be set equal to the internal 
momentum or other scales such as k'^ or A^. The constraints for the exponents 7^, 
(i arise from RG rescaling invariance. 

4.4 Numerical results 

4.4.1 Perturbative initial conditions and bare vertices 

Here we solve the flow from the perturbative UV towards the IR, imposing one-loop im- 
proved initial conditions. For illustration, we start the flow at A = Mz = 91.187 GeV with 
a{Mz) = 0.118. 

Our results for the bare-vertex truncation are shown in Fig. |21 In agreement with our 
analysis of Sect. 14. 2[ both dressing functions increase towards the IR. After a few orders 
of magnitude of perturbative running, the flow in fact becomes singular at a finite scale 
fcsing where the gluon dressing diverges.^ As a consequence, the coupling a also diverges 
reminiscent of the Landau-pole singularity of perturbation theory. 

^This gluon divergence should not be confused with a behavior of the gluon propagator, as conjec- 
tured earlier in the context of infrared slavery. The divergence of the flow at fcsing implies that the gluon 
propagator diverges at a finite value of p^. Hence this divergence signals a breakdown of the truncation 
and not the onset of a behavior. 



V,A = G^^ip')G^'ir')Z'^-ip')Z'^-ir'), 71 + 72 = 2, Ci + C2 



-1 



(35) 
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Figure 2: Numerical results for ghost and gluon dressing functions and the running cou- 
pling employing the bare vertex truncation. Left panel: bare 3-gluon vertex (a). Right 
panel: modified bare 3-gluon vertex (b) with p decreased from 1 to 0.2. Both flows become 
singular at a finite value of /c, signaling the breakdown of the oversimplified truncation. 
The coupling exhibits a singularity of Landau-pole type. Other multiplicative positive 
3-gluon vertex dressings lead to similar results. 



This singular behavior of the gluon is mainly driven by the gluon loop involving the 
3-gluon vertex. Since the latter is expected to be suppressed in the IR, let us use the 
modified bare 3-gluon vertex with a suppression parameter of p < 1. On the left panel of 
Fig. 121 we show the flow for this truncation with p gradually decreased from 1 to 0.2 for 
decreasing k. Though the gluon dressing now remains finite, it is the ghost dressing that 
diverges at finite k and thereby triggers a singularity of Landau-pole type. 

These observations are in agreement with the analysis of Sect. 14.21 and stress the fact 
that a truncation with a multiplicative positive dressing of the 3-gluon vertex is not capable 
of describing the flow from the perturbative UV to the IR power-law asymptotics. 

4.4.2 Perturbative initial conditions with dressed 4-gluon vertex 

Let us now include a momentum-dependent 4- vertex dressing of the type (e) in Eq. (jH^ . 
In agreement with our observations of Sect. 14.21 the truncation now has the potential to 
describe an IR suppression of the gluon towards the expected IR asymptotics. 

In Fig. El we display a typical solution for the gluon and ghost dressings, exhibiting 
IR ghost enhancement and gluon suppression. All solutions look qualitatively similar for 
vertex parameters 7i,Ci > in Eq. (j35|) . We also observe a quantitative stability for 
7i = Oil) and Ci = 0.1 . . .0.5. (Fig. His computed with 71 = 1.1, Ci = 1/2, = and 
a modified bare 3-gluon vertex (b) with p decreasing from 1 to 1/2 in the mid-momentum 
regime in order to avoid a Landau-type singularity). 

Nevertheless, the vertex ansatz fl35j) appears to be too simple for establishing a full 



18 



o.s' , , I 

10"^ 10'' 10° lo' 10" 10' 10* 10^ 

p' [Gev'] 

Figure 3: Numerical results for the ghost and gluon dressing functions, including a dressed 
4-gluon vertex (e). This truncation is capable of describing ghost enhancement and gluon 
suppression in the IR. (Vertex parameters: 71 = 1.1, (i = 1/2, = A^, and p = 1 . . . 1/2.) 

UV-to-IR connection. Although Landau-type singularities can be avoided over the full 
momentum window that we have studied, we have not discovered an IR asymptotics with 
a high degree of universality, such as a clear signal of stable power laws for the dressings. 

4.4.3 IR flow analysis 

Let us now concentrate on the flow of the propagators towards the IR asymptotics. DSE as 
well as RG flow equation studies have demonstrated that the power-law behavior discussed 
in the introduction is a self-consistent solution of the dynamical equations in the IR, 
implying an IR stable fixed point for the gauge coupling. The following study is devoted 
to an investigation of the domain of attractivity of this fixed point. For this, we assert 
that the fluctuations from the UV down to an intermediate scale A ~ IGeV have already 
modified the dressing functions as compared to their perturbative form in an a priori 
unknown manner.^ By assuming various initial conditions for the dressing functions at 
intermediate k scales, we can check whether the flow connects a particular mid-scale initial 
condition with the IR fixed point regime and power laws. In this way, we can analyze the 
mid-momentum requirements that facilitate a full UV - IR connection. This particularly 
provides information about the physical mechanisms that have to be triggered by the 
(yet unknown) full vertices. Throughout this subsection, we use bare ghost and 4-gluon 
vertices - a truncation that is sufficient for the expected IR asymptotics, owing to ghost 
dominance in the gluon equation. Beyond this, we vary the 3-gluon vertex in order to 
study its influence on the approach to the IR. 

^In this analysis, we choose the scale by matching the perturbative one-loop expression for a{p^) to our 
results and employ the experimental input a{Mz) — 0.118 (even though we neglect dynamical quarks). 
The intrinsic uncertainty of such a procedure is large, but irrelevant for our purposes here. 
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[GeV^] [GeV^] 

Figure 4: Numerical results for the ghost and gluon dressing functions in the minimally 
dressed vertex truncation (c), Eq. (jH^ . employing a gluon dressing function with a scale, 
(|3ip. and a constant ghost as boundary conditions at an intermediate momentum scale. 
For an initial gluon power in a small but finite interval around r^^ ~ 1, the ghost shows a 
clear signal for IR power-law asymptotics with /t ~ 0.52. The right panel magnifies the IR 
asymptotics of the ghost of the left panel, demonstrating the failure of the ansdtze with 
TjA = 0.8 or TjA = 1.2, which are outside the attractive region of the IR-power law solution. 
(The vertical dotted line indicates the lower end /cir, of the window of k integration.) 

Two initial conditions for the ghost and gluon dressing functions at = A reflecting 
opposite situations in the mid- momentum regime have been given in Eq. ()30|) and (jHT|) . In 
addition to being related to A, both boundary conditions contain a scale L appearing either 
in the initial ghost or gluon dressing function. Certainly, Yang- Mills theory is governed by 
only one scale, which is Aqcd- The relation between A and Aqcd is uniquely determined by 
the RG, once the coupling g\ is fixed at A. Such a unique relation also has to exist for L, i.e., 
L = const. X Aqcd- To guarantee that our initial condition represents a valid approximation 
to Yang-Mills theory at our initial scale A, we have to determine this constant. Otherwise 
our system describes a different theory with two independent scales. 

For solving this problem, we note that, once we have found the right value of L, the 
flow in the deep IR will not develop yet another scale and the dressings will depend only 
on p^/L^ ~ p^/Aqcd- On the other hand, a separate dependence of our solutions on L 
and A indicate the failure of the initial condition. In practice, instead of tuning L/A for 
fixed qa, we keep L/A fixed, say L/A = 0.01, and fine-tune the value of the initial coupling 
g\ such that no further scale arises in our solutions in the deep IR. Examples of such an 
additional but unphysical scale are rapid changes, steps or singularities in the dressings. 

Let us first investigate the initial condition (j30|) . which assumes that the ghost dressing 
has already developed a scale, whereas the gluon is taken as perturbative for all momenta 
p'^. With these boundary conditions we find that the flow runs into a singularity at finite 
ksing independently of the initial ghost exponent rjc or the value of L. The mechanism 
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Figure 5: Numerical results for the ghost and gluon dressing functions, employing a 
gluon dressing function with a scale, (jSH) and a constant ghost as boundary conditions at 
an intermediate momentum scale. The left panel shows the result for the modified bare 
vertex truncation (b) with p = 0.5, a "ghost-loop-only" truncation (p = 0) is displayed on 
the right panel. Both truncations lead to IR power law asymptotics with k ~ 0.52 if the 
initial gluon power is in a small but finite interval around rjA — ^■ 



responsible for such a behavior has been analysed in Sect. 14.21 in the truncation considered 
here, the gluon cannot develop an IR suppression with — > 0. Thus an unsuppressed gluon 
at a mid- momentum scale k = A drives the fiow away from the IR fixed point. 

This is different for the initial condition (jHT|) . where we assume that the gluon has 
already developed a suppression with scale L at k = A and the ghost remained perturbative. 
In Fig. m we plot the solution of the flow equation for the initial conditions with initial 
gluon exponent 7]j^ = 0.8,1.0,1.2 for the minimally dressed 3-gluon vertex (c), Eq. (j34j) . 
The regulator scale k has been integrated down to a value /cir that is indicated by the 
vertical dotted line on the left-hand side of the plots. For all three values of rjA, the ghost 
dressings are strongly enhanced in the IR, whereas the gluon develops a more enhanced 
peak around L. 

For rjA = ^ (solid lines), the ghost dressing exhibits a clear power law, G{p'^) ~ {jp')~'^, 
for all momentum values kj^ ^ ^ L^- We identify a ghost exponent of 

K ~ 0.52. (36) 

The gluon develops a power slightly stronger than its initial power rjA- This behavior is 
universal for initial powers rjA in a small but flnite interval around r/^ — 1, i-e., for boundary 
conditions in this interval the IR power k does not depend on rjA- As will be discussed 
in more detail below, this result corresponds to the self-consistent solutions found, for 
instance, in the IR analysis of DSEs. 

For 7] = 1.2 (dot-dashed lines), a power law in the ghost dressing develops, but not over 
the full k integration range. At some scale k"^ ^ /cjj^, the ghost gets stuck and becomes 
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essentially perturbative again. Apparently, the gluon is too strongly suppressed to control 
the build up of a power law to arbitrarily small momenta. Since this solution exhibits a 
second independent scale where the ghost becomes perturbative again, we conclude that 
this ansatz with tja = 1-2 fails to describe the IR of Yang-Mills theory. 

For TjA = 0.8, the ghost first builds up a power law, but then deviates in the IR and runs 
into a singularity (the plot shows a curve with g/^ chosen such as to allow for a maximal IR 
extension). For this set of initial conditions, the gluon is not suppressed strongly enough 
in the IR to prevent the same singularity as observed frequently above. We conclude that 
also this initial condition with rjA = 0.8 fails to represent the physics of the IR sector. 

These features remain qualitatively the same for the various 3-gluon vertices, unless 
this vertex is not too strong at intermediate scales. In Fig. we show results for the 
modified bare vertex (b) with p = 0.5 and a "ghost-loop-only" truncation (p = 0). In 
agreement with the literature, this reveals the ghost loop as the dominant IR structure, in 
support of the Kugo-Ojima confinement scenario. 

5 Summary and conclusions 

We have performed a study of the vertex expansion of the quantum effective action of 
Landau-gauge Yang-Mills theory in the framework of the exact renormalization group. We 
have truncated the vertex expansion at the lowest nontrivial order and concentrated on 
the behavior of the ghost and gluon propagators. This level of truncation already contains 
important information about large-distance physics along the lines of the Kugo-Ojima 
confinement criterion. 

We identify three different momentum regimes in our investigation: a perturbative 
high-momentum region, a nonperturbative regime around ~ (9(l)GeV^ presumably 
dominated by gauge field fluctuations and a regime in the deep IR dominated by the ghost 
degrees of freedom. Simple truncations of the RG flow equation cannot only easily deal 
with the perturbative regime, but give conclusive answers for the deep IR. In the literature, 
a power-law behavior of the dressings has been self-consistently determined in the deep IR, 
which is related to an IR fixed point of the running coupling and realizes the Kugo-Ojima 
confinement criterion [HI 13 IHl IHO] • 

In addition to these studies, our work analyzes the flow towards this fixed point and 
explores its domain of attractivity. For this, we employ a set of mid-scale initial condi- 
tions for the propagator dressings that may arise from the flow, once the fluctuations at 
intermediate momenta have also been properly integrated out. Within this set, our results 
demonstrate that a mass-like structure in the gluon dressing (characterized by an initial 
gluon power in a small but finite interval around r]A — 1 and a scale L) is a necessary 
prerequisite for approaching the power-law asymptotics in the IR. Other tested mid-scale 
initial conditions lead to either a singular or a trivial behavior of the propagator dressings. 
For the mass-like structure, the resulting flow shows a universal behavior, being largely 
independent of the initial details, such as further properties of the initial conditions or the 
form of the gluonic vertices. 
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As a result, we observe a ghost exponent of k ~ 0.52. Even though the corresponding 
value of the DSE IR analysis is k ^ 0.595, these results are in satisfactory agreement, 
since the exponents are regulator dependent in the present truncation. This regulator 
dependence can quantitatively be studied in a self-consistent IR analysis in the flow equa- 
tion framework revealing k ~ 0.539 for the regulator used in the present work |3H]; 
the difference to our result obtained by solving the flow equation gives a measure for the 
numerical accuracy of our procedure. 

We furthermore have searched for global solutions connecting both asymptotic ends 
of the momentum range. Here the mid-momentum regime appears most difficult and 
represents a delicate problem. We show that oversimplified truncations based on bare 
vertices lead to Landau-type singularities in the flow, signaling the importance of higher- 
order correlations. Dynamical information encoded in dressed (or "running") vertices is 
required in this mid-momentum regime. 

In the present work, we have modeled this dynamical information by supplementing 
the bare-vertex index structures with (positive) momentum-dependent dressing functions. 
Within these limits, we observe that a truncation on the 3-point level typically leads to 
ghost and gluon enhancement in the IR, and is therefore insufficient to describe the IR 
suppression of the gluon dressing function, as, e.g., observed on the lattice [121 CH IT3] . 
Once suitable momentum-dependent 4-point dressings are taken into account, the gluon 
does indeed exhibit the expected IR suppression. We have demonstrated this property by 
garnishing the 4-gluon vertex with generic powers of the propagator dressing functions. 
However, even though the resulting flow remains stable over all momentum scales of the 
numerical integration, the resulting propagators in the deep IR depend strongly on the 
details of the modeling with no sign of the expected universal behavior. We conclude that 
in the truncations considered, important dynamical information from the mid-momentum 
regime is still missing. 

As a main conclusion, we interpret our results as follows: the IR power-law asymp- 
totics of the propagators is triggered (though not generated!) by a particular dynamics at 
intermediate momenta. If this dynamics gives rise to a mass-like structure in the gluon 
propagator, the Kugo-Ojima branch acts as an IR attractive fixed point of the flow, re- 
sulting in ghost enhancement and gluon suppression in a universal manner. Our results 
illuminate the role played by the mid- momentum regime, but, apart from some qualitative 
insight, the physical mechanisms potentially governed by higher-order correlations are still 
poorly understood. 

Concerning quantitative estimates made in this paper, let us note that we have ex- 
pressed all momenta in physical units (GeV) by fixing the flow at the Z mass, Mz- We 
have done so mainly for illustrative purposes. We stress that resulting scales in the nonper- 
turbative sector can sizably be modified by changes in the vertex dressings and, of course, 
upon the inclusion of dynamical quarks. 

From a conceptual perspective, an important ingredient for our findings is the modified 

^"^In fact, since the flow has to build up the full ghost power from an initially flat dressing, numerical 
errors generically accumulate such as to reduce the final value of n. 
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Ward-Takahashi identity (mWTI). Beyond the dynamical information contained in the 
flow equation, the mWTI acts as an additional constraint that encodes gauge invariance 
even in the presence of a gauge-symmetry breaking regulator. In the present case, we 
have employed the mWTI for controlling the running of the gluon mass rrik at finite k; 
this turns the gluon mass from a seemingly independent and relevant RG parameter into a 
dependent and RG irrelevant quantity. In the language of loop integrations, the gluon mass 
corresponds to quadratic divergencies arising in non-gauge-invariant regularizations for 
which appropriate subtraction procedures have to be defined. We have verified that naive 
subtraction procedures not obeying the mWTI generically generate an artificial mass scale 
in the flow. This in turn also leads to a global solution which smoothly interconnects the 
perturbative and the IR power-law behavior. However it does so for the wrong mechanism. 
The scale L discussed above certainly should be generated dynamically and not by the 
regularization procedure. 

Future work should be devoted to an understanding of the dynamical mechanisms 
that initiate gluon suppression, as required for approaching the IR attractive domain of 
the Kugo-Ojima branch. A possible route is the inclusion of the next order in the vertex 
expansion, pursuing a dynamical calculation of vertices. Promising work in this direction in 
the flow equation framework has been done in ^\ , where an attempt was made to compute 
momentum-dependent 4-point vertices. The resulting quantitative deviations from other 
results in the literature may be associated with the use of an angular approximation and the 
neglect of the 3-gluon vertex, but improvements are certainly straightforwardly amenable. 

Finally, it remains to be understood how the vertex expansion can be connected with 
another nonperturbative truncation scheme: the expansion of the effective action in field- 
strength invariants [21] • Whereas the vertex expansion concentrates on the flow of opera- 
tors with a small number of fields but full momentum dependencies, the expansion in terms 
of field-strength invariants can deal with an infinite number of fields but usually neglects 
momentum dependencies. It is reassuring that an IR fixed point in the coupling has also 
been observed in the latter expansion of the effective action ^B] ? in particular, since both 
approaches exploit a non-renormalization theorem in order to define the running coupling. 
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A Integral kernels 

Here we give the explicit formulae for the momentum integral kernels of the propagator 
equations (fT7|) and (fTHj). In the Landau gauge, only the kernels involving transversal gluons 
are of primary importance. For the bare-vertex truncation, they read 



Q^'^(x,y,n) = -{l-u") 
2{D-X 



x + y -\ 1 



xy ( \ — V? 



z \ D-1 



(A.l) 



= (A.2) 



' 1 — ti^) X y 



Q^'^(x,y,n) = ^ j^^, (A.3) 

where the variable z on the right-hand side has to be understood as an abbreviation of 
z = X + y — lu^fxy. 

B Regulator properties 

For the regulator shape function specified in Eq. (jD 



r(x) = rAc(:r) = (B.l) 

X 

the following identities hold: 

xr{x) = {\-x)B{\-x), -x^r'{x) = d{\- X). (B.2) 
These are required in Eq. (fT^ . The regularized propagators finally read 

1 _ Q{\-x) Q{x-\) 
Pa(x) 1 + m^//c^ x -|- m^/Zc^ ' 

+ — -. (B.3) 



VcKx) X 

C RG rescaling 

Let us study the system under (finite) RG rescalings of the renormalized fields of the form 

^ Zl^l^A^, Z-JI^ (C.l) 

where Z^^ Zc account for a possibly fc-dependent RG scaling. Invariance of the action 
under this scaling transformation implies that the n-point vertices scale as 

p(n+m) p(n+m) /p r)\ 

^ A C A^ ' V / 
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The flow equation is invariant for regulators that scale as 

Ra 2, A Rai Rc Rc- (C-3) 

As a special case of Eq. ()(I2|1 . let us note that the propagator dressing functions scale as 

G{x) G{x), Z{x) Z[x). (C.4) 

Owing to its nonperturbative definition the renormalized gauge coupling scales as 

g^-.ZlZj.g^ (C.5) 

such that g^G^Z is invariant as required. In Eq. (|lUj). we defined the vertex dressings V 
by scaling out the coupling, T^^^ = gV^^\ T^^^ = g'^V^^^] hence the vertex dressings have 
to scale as 

^CAC ~^ ^CACi ^3A — ^ ZaZ^^ Vsa, V^a ^ ZaZ'^^ V^a- (C-6) 

These properties tell us, for instance, that the naive bare vertex truncation Vsa 1 
is generally not compatible with RG scaling, and therefore only correct if the fields are 
normalized canonically. 

D Numerical procedure 

For the numerical computations, we have to specify the regulator shape function introduced 
in Eq. (fT^ . Here we use the one proposed in |35|, 

rA,c{£) = lil-x)9il-x). (D.7) 

The step function in Eq. ()D.7|1 leads to a more localized support of the loop momentum 
integrations and thereby simplifies the required algorithms. 

For the actual numerical integration of the flow from the starting scale A to an IR scale 
kin, we employ a fifth-order Runge-Kutta formula with an embedded fourth-order formula 
to estimate integration errors. We have also tested Euler's method and found it neither 
accurate nor stable enough for the present problem. The loop integrals on the right- 
hand side of the flow equations are carried out numerically; no angular approximations 
are made. To this end, the dressing functions have to be represented over the whole 
momentum range. We choose a Chebyshev expansion on a logarithmic momentum grid, 

G [e^, Ap] with ep -C /c|^ and ^ A. An expansion of the dressing functions in terms 
of the first Chebyshev polynomials is exact for those A^ momenta that are zeros of 
the (A^ -|- l)st polynomial. Thus, taking these zeros as external momenta p'^ in the RG 
equations ensures maximal numerical stability. A typically chosen value for A^ is A^ = 45. 

^^Whereas the nonanalyticity of the 9 function can become problematic in high-order derivative expan- 
sions jSni, there is no conflict with the vertex expansion and the properties of this shape function wiU 
remain numerically useful also to higher orders. 
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For momenta far below e^, we extrapolate the Chebyshev expansion employing either a 
power-law fit or a smooth continuation to a constant value; the latter is justified by the 
presence of the regulator. We have confirmed that our results do not depend on the details 
of this extrapolation. In general, the advantage of the Chebyshev representation over the 
parametrization method used in jSEI is that one does not rely on a specific form of fit 
functions in the interval [e^, A^] which could possibly bias the results. 
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